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Abstract. - The role of memory is crucial in determining the properties of many dynamical pro- 
cesses in statistical physics. We show that the simple addition of memory, in the form of noise 
reduction, modifies the overall scaling behavior of the voter model, introducing an effective surface 
tension analogous to that recently observed in memory-based models of social dynamics. The nu- 
merical results for low-dimensional lattices show a scaling behavior in good agreement with usual 
Cahn-Allen curvature-driven coarsening, even though slower preasymptotic regimes may be ob- 
served depending on the memory properties. Simple arguments and a mean-field analysis provide 
an explanation for the observed behavior that clarifies the origin of surface tension and the mech- 
anism underlying the coarsening process. 



The theory of phase-ordering kinetics [1] describes the scaling behavior of domain coars- 
ening phenomena occurring in non-equilibrium statistical physics, with applications ranging 
from social sciences [2] to ecology [3] . The way in which a system orders starting from a dis- 
ordered phase follows different paths depending on generic properties (e.g. symmetries, the 
type of interactions, etc.) that allow to identify specific universality classes of nonequilib- 
rium phase-ordering dynamics. In particular, the important family of 'kinetic Ising models' 
with ^-symmetry and no strict conservation of the order parameter displays two different 
ordering behaviors, depending on whether a surface tension exists or not. If this is the case 
(as in Glauber dynamics, for example), coarsening is curvature-driven [4]: the domain length 
scale grows following a power law, £(t) ~ t 1 / 2 in any dimension [1]. Other models, such as 
the voter model [5] , exhibit domain coarsening without surface tension: dynamics is driven 
by interfacial noise [4]. In this case £(t) grows as t 1 / 2 in one dimension, logarithmically in 
d = 2 and it does not diverge for d > 2 [6] . It is also possible to define a family of kinetic 
Ising models interpolating between zero-temperature Glauber dynamics (TOGD) and the 
voter model (VM) [4,7,8]: in these models, the transition rates depend on two parameters 
measuring the strength of bulk and interfacial noise. When bulk noise is present the coars- 
ening is curvature-driven. When only interfacial noise exists, coarsening is voter-like for a 
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broad class of models, cither characterized by global conservation of magnetization or by 
nonconserving Z 2 -symmetric dynamical rules [4]. 

In this Letter, wc introduce a simple model, consisting of a voter-like dynamics modified 
by a noise-reduction procedure. The main motivation of this work consists in understanding 
the mechanism governing the ordering process of some recently proposed models of social dy- 
namics, such as the Naming Game [9,10] and the Wang-Minett model of bilingualism [11,12]. 
In both these models, the pairwise interaction rules are reminiscent of the voter dynamics, 
but the update follows a two-steps process involving memory effects. Surprinsingly, on 
d-dimcnsional lattices, these models undergo a domain coarsening dynamics following ap- 
proximately the same temporal laws of curvature-driven models with non-conserved order 
parameter. In order to elucidate the origin of this effective curvature-driven coarsening, 
we define a minimal model, that we call "noise-reduced voter model" (NRVM), in which a 
noise-reduction prescription is added to the usual interaction rule of the VM, introducing a 
memory effect. As it will be clear below, the introduction of a noise-reduction prescription 
does not give rise to bulk noise in the spin dynamics. Instead, it turns out that the addi- 
tional ingredient gives rise to an effective, memory-induced, surface tension, that generates 
in low dimension a curvature-driven coarsening process. 

Noise reduction has been used in the study of surface growth [13] and diffusion-limited 
aggregation [14] as a method to speed up the approach to asymptotic scaling. It amounts 
to attaching to each site of the system a counter recording how many times a particle has 
touched the site. Only when the counter reaches a given threshold value the site becomes 
active and the particle sticks to it. We introduce the noise reduction in the framework of 
the VM in the following way: each site i is endowed with a spin variable, Sj, as well as with 
a pair of counters Cf and C[~ . At each time step, a site i and one of its neighbors j arc 
selected at random to interact. The counter associated with the spin value of the selected 
neighbor is increased by one: C^ 3 — > C* 3 + 1. At odds with the normal VM dynamics, the 
spin variable s, is not necessarily modified upon interaction. It is updated only when one of 
the counters attains an a priori fixed threshold value r (e.g. if the positive counter reaches 
the value C 4 + = r then s, is set equal to +1). In this way the interaction between spins is 
indirect, being mediated by the dynamics of counters. Two possible update rules for the 
counters are considered: when the spin at site i flips, only the counter that has reached r is 
set to zero (model 1), or instead both counters are reset (model 2). 

For r = 1 both models coincide with the usual voter dynamics, whose behavior is well 
known in any dimension [6] and is characterized by the absence of surface-tension [4] . For r > 
1, one expects an effective surface-tension to arise: the repeated interaction with neighbors, 
needed before a spin update is actually performed, implies that the counters sample the 
local field of each site. According to this argument, in the long run, local field imbalances 
are accumulated, pushing the system evolution rule towards the TOGD. In fact, the NRVM 
does not simply interpolate between the VM and the TOGD: the transition rates encode 
the memory on how the spin pattern has been reached, thus it is easy to find examples of 
configurations such that the rates are not intermediate between regular voter and Glauber 
dynamics. 

In order to ground more firmly the intuition about the existence of a surface tension, we 
have checked the presence of surface tension in a two-dimensional lattice model by studying 
the evolution of a circular droplet in a large system of N = L 2 sites. Figure 1 shows that, 
as for systems with curvature-driven coarsening [1], the volume of the droplet decreases 
linearly with time, i.e. the typical size of a domain decreases as \ft. The figure refers to 
model 1 dynamics, but a similar behavior is observed even for model 2. Moreover, rescaling 
time with r the curves for different values of r are very close, revealing that changing the 
threshold value produces only a linear temporal rescaling in the ordering process. 

An explanation of the linear rescaling can be argued from the analysis of the equation 
for the dynamics of counters. The probability that a spin Si flips is given by the joint 
probability that a neighboring site of opposite spin is chosen and the counter C^ Si = r — 1. 
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Fig. 1: Illustration of the presence of surface tension in the two-dimensional voter model with noise 
reduction. The initial condition is a configuration (L — 500) with a droplet of radius Ro = 200 
surrounded by a sea of opposite spins. We have monitored the normalized area nR 2 {t)/L 2 as 
a function of time for model 1. This quantity decreases linearly for r > 1, as expected in the 
presence of surface-tension [1]. Note that the linear trend is clear even if the data refer to a single 
realization. Time is rescaled by r to show that the noise-reduction parameter simply rescales the 
temporal evolution. In the inset, the average time r needed to dissolve the droplet is plotted as a 
function of r for both models: r scales linearly for model 1 and with an effective exponent 1.1 for 
model 2. 



On a (i-dimensional lattice, the probability to choose a neighbor of i with opposite spin value 
is W(si) = | 1— ^ s iSjeV(i) s i ' wnere is the neighborhood of i. In model 1, the 
dynamics of the two counters at a site i are decoupled, thus the probability p^ Si (t) that the 
counter C 2 ±Sl is equal to a value n at time t is defined by the equation, 

§- t Pi s ^t) = - w( TSl ) , (i) 

with periodic boundary conditions on the variable n. The solution of this equation requires 
the knowledge of the probability W(si) at all times, i.e. of the whole dynamics but, formally 
solving the equation in the Fourier space, we can write 

Pr-M = - I 1 + Y, cos ( 2lTk / r + sin(27rfc/r)Vtf (*)) e -[i-«»P**/r)]Wf (t) j ; (2 ) 



k=l 



where Wf{t) = /„ W(±Si(t'))dt' . Note that the sum in the r.h.s. of Eq. 2 tends to when 
the integral is much larger than r, while it tends to -1 when it is much smaller. For a spin in 
the bulk of a domain W(—si) = 1, while W(si) — 0. Hence, in the former case the integral 
W~(t) diverges and p^i^t) — > 1/r, while in the latter the integral (t) is finite, so that 
Pr-i(t) — > 0. As a consequence the effect of varying r ^> 1 can be reabsorbed in a temporal 
rescaling t — * t/r (see the inset of Fig. 1). In order to extend the analysis to model 2, 
one should consider the equation for the joint probability 9^ i n _ (t) that the counters Cf Si 
assume some values n + and n_, since the counters are now coupled. We will consider such 
approach later for the simpler case of a mean-field model with r = 2. However, we find 
numerically that in d = 2 the scaling is not far from the one for model 1. In d = 3, changing 
r in simulations apparently gives more than a simple rescaling of time (see below), but this 
is probably just a preasymptotic effect. 

Another remark arises considering the master equation for the probability distribution 
P(S,t) of having a certain d-dimensional configuration S = {si} at the time t in model 1, 

j t P{s,t) = ]T [wi-sjp^txuit) - w( Si )p(s,t)p;f_\(t)] , (3) 
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Fig. 2: Temporal evolution of the density of interfaces nA(t) in the voter model with noise reduction 
on a two-dimensional lattice of size L = 5000 (left) and on a three-dimensional lattice of size L — 300 
(right). The straight lines, reported for reference, have the slope 1/2 expected for curvature-driven 
coarsening. 



where P(S l ,t) is the probability of the configuration S with the only spin Si flipped. A 
similar equation holds for model 2. From Eq. 3, we can write the evolution equation for the 
average quantities, such as the local magnetization (sj) 

^ <x (p^tjVfr) , (4) 

in which p~l\{t) evolves according to Eq. 1. Eq. 4 is similar to the analogue equation for 
the usual voter model [6], but the presence of the counter probability p~ Si breaks down the 
conservation of the average magnetization, which is one of the crucial features of the voter 
model on regular lattices. 

The presence of an effective surface tension is expected to induce a curvature-driven 
coarsening, resulting in a phase-ordering process with scaling properties belonging to the 
Cahn- Allen universality class [1]. This has been checked directly by studying numerically 
the coarsening dynamics in d = 2 and d — 3 starting from a fully disordered random 
configuration with zero magnetization. The temporal evolution of the density of interfaces, 
riA^t), reported in Fig. 2, shows in general a power-law decay oc t~ x l z for all cases. 

For model 1, after the trivial initial transient, the effective exponent sets to about 0.45 in 
d = 2 and close to the value 1/2 typical of the Cahn- Allen universality class in d = 3. 
Remarkably, the value 0.45 coincides with the value found recently by Castellb ct al. [12] 
in another modified voter model. We cannot ascertain numerically whether the discrepancy 
with 1/2 is the consequence of logarithmic corrections or a truly asymptotic exponent. For 
model 2, the violation of the curvature-driven decay (exponent 1/2) is much stronger. In 
d = 3 a numerical fit gives an effective exponent that decreases with r: 0.39 for r = 20 and 
0.33 for r = 100. The origin of this change can be understood by visual inspection of the 
pattern of coarsening domains (Fig. 3). For large r the interfaces between the two phases 
tend to be faceted up to a length £ c (r), that grows as r is increased. At the beginning of the 
ordering process, when the typical size of domains £(t) is smaller than £ c (r), the coarsening 
process is strongly perturbed (and slowed down) by the faceted nature of interfaces. This 
explains the reduced value of the effective exponent. Later on, when the scale of domains 
becomes sufficiently large (£(t) £c(f)) the short scale structure of domain boundaries 
does not matter anymore, interfaces will be rough and therefore we expect curvature-driven 
coarsening (with its associated exponent 1/2) to rule. This picture applies well to the two- 
dimensional case. The effective exponents found for large r in d = 3 are probably due to an 
extremely long crossover to the asymptotic regime, although present numerical evidence is 
not clear-cut in this respect. 

To understand the dynamical interplay between the voter dynamics and the nonequilib- 
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Fig. 3: Snapshot of a spin configuration for model 2, r — 2 (left) and r = 1000 (right) showing 
that domain interfaces tend to be faceted for large r. The two configurations have comparable total 
length of the interface. 



rium surface-tension, it is useful to perform a mean-field analysis, considering the NRVM 
on the complete graph. Following Ref. [15], we derive a master equation for the probability 
P(m, t) that the magnetization takes the value m at time t. In a mean-field framework, 
single-site counters are well represented by the average values over the whole system. 
We call pfi (t) the mean- field probabilities that the average counter value is equal to n 
at time t. In a temporal step, the magnetization can change of a value ±2/N or remain 
constant. The corresponding probabilities appearing in the master equation are, 

| 2 1 1 + m 1 — m ± 
Pvoh s Tfi > Tfl zb -j^ ? — — —P r -i(t). (5) 

Writing down the master equation, using Eq. (5), expanding the expressions up to the second 
order in l/N (with N — > oo), and neglecting higher order terms, we get 

^ P K*)= 2^[ftl(*)-Pr-l(*)]^[(l-"* 2 )^K*)] 

+ ^2 [Pr-l(*) +JV-l(*)] ^2 [(1 ~»i 2 )P(m,i)] . (6) 

Two other equations for the counter probabilities Pn(t) are required to close the equation 
for P(m,t). For model 1, the equations have the same form of Eq. 1, while in the case of 
model 2 one should start from the evolution equation for the mean-field joint probability 
<?n + ,n_ (t), noting that p r _ x (t) — ^ n± <?n + ,n_ (t). The solution of the coupled set of equations 
is, in general, beyond reach, but it is instructive to look at Eq. 6: apart from the counters 
dependent prefactors, the second term is the one present in the equation for the VM [15], 
while the first is instead the one found for TOGD [16], expressing the existence of a drift, 
due to the imbalance between the counters. The interplay of the two terms reflects the 
competition between curvature-driven and noise-driven coarsening. 

The rate equations for the mean-field magnetization m(t) can be similarly derived; let 
us elucidate the simple cases for r = 2. In model 1, the following system of rate equations 
holds 

m(t) = 1 ~ n f> 1 [Pt(t)-pz(t)] , (7) 
p-5 
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Fig. 4: Temporal evolution of the magnetization m(t) for model 1 (left) and model 2 (right) on 
a complete graph, using rate equations. Two curves for different initial magnetization mo reach 
different stationary states in model f . 

P f(t) = '-^[i-iptm . (8) 

The case of model 2 is more complicate, since we need to compute the four coupled evolution 
equations for q n+ ,n-(t). The corresponding system of rate equations reads 

1 _ ™(+\ 2 

(9) 
(10) 

(11) 
(12) 

(13) 



m(t) = 


1 - m(tf 
2 


kl,o{t) - 90,l(*)] , 




9o,o(*) = 


-?o,o(*) + 


1 + m(t) , . 1 
2-^<?i,o(t) + - 




9i,o(*) = 


-9i,o(*) + 


1 + m(t) ^ 
2 9o,o (*) , 




9o,i(<) = 


-<fo,i(*) + 


1 - m(t) 

2-^«0,o(*) , 




(<) = 


-qx,i(t) + 


1 + m(t) , . 1 
2 A2 «),i(*) + - 


— m(t) , . 
- T Ai 9i,o(«) , 



Linear stability analysis indicates that the two models present different stationary solu- 
tions. For model 2, there are only two possible stable stationary solutions, {m = — 1, go,o = 
1/2,51,0 = 0,go,i = l/2,gi,i = 0}, and {m = +l,£?o,o = V 2 >9i,o = l/2,g ,i = 0, = 0}, 
i.e. the system completely orders into a ferromagnetic phase. The solution with m = is 
instead unstable. Hence, in the model with double counter update, the mean-field dynamics 
is completely governed by the drift term. On the contrary, in the model 1, fluctuations 
play a relevant role. The possible stationary states are p± = 1/2 for all possible values of 
m, and m = 1 and = 1/2 for any value of pj" < 1/2 (or the symmetric case m = — 1, 
Pi = 1/2, pi < 1/2). The existence of multiple stable stationary states points to the ex- 
istence of a transition depending on the initial value of the magnetization to(0). Indeed, 
solving numerically Eqs. (7-8) we find 

f ±m(0)/m* for |m(0)| < m* , . 

~\ ±1 for |m(0)| > m* [ > 

while the difference pf — p{ is below the transition and finite for |m(0)| > m* (Fig. 4). 
The transition point is m* « 0.6382 . . .. 

If one performs simulations on the complete graph, the sharp transition is actually smeared 
out by fluctuations. Nevertheless we can observe a rather abrupt crossover between a regime 
where pf — Pi is exponentially small (for |m(0)| < 0.22 . . .) to a regime with finite pf — p± , 
whereas asymptotically \m\ is equal to 1 in both cases. The shift in the position of this 
transition is due to the diffusive contribution of the second term in l.h.s. of Eq. 6, that 
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is neglected in the rate equations. For small |m(0)| the fully ordered state is not reached 
because the initial condition is not sufficiently asymmetric and the drift, which decreases 
during the evolution (both and tend asymptotically to 1/2), is not strong enough to 
fully order the system. If |m(0)| is large enough instead, the counter imbalance is sufficient 
to drive the system to full order, keeping an asymptotic nonvanishing drift. For r > 2 the 
general phenomenology remains similar, with the difference that the counters probabilities 
may exhibit oscillations that in turn induce a nonmonotonic behavior of the magnetization. 

In conclusion, we have introduced a noise-reduced voter model, in which memory effects 
induce an effective surface tension. While in finite dimension this implies that the coarsening 
process is always curvature-driven, on a fully connected graph it is possible to observe a 
crossover between voter-like and zero temperature Glauber behavior. We have provided an 
explanation of the mechanism generating surface tension, and numerical results supporting 
the thesis that the NRVM ordering process belongs to the Cahn- Allen universality class of 
coarsening with non-conserved order parameter [1]. Further investigations on model 2 are 
required in order to clarify the asymptotic value of the coarsening exponent for large r. The 
present analysis is consistent with recent results obtained by Castello et al. [12], showing 
that the existence of an intermediate state (here due to the introduction of memory) is 
relevant in modifying the scaling properties of the VM. In order to give further insight on 
more realistic situations, such as those involving memory-driven non-equilibrium ordering 
dynamics on social networks, it would be interesting to study the behavior of the model on 
systems with heterogeneous topologies. We expect that the competition of the two types of 
kinetic mechanisms, i.e. curvature-driven coarsening and intcrfacial noise, may give rise to 
nontrivial ordering phenomena. 
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